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The concept of the electron locaUzation function (ELF) is extended to two-dimensional (2D) 
electron systems. We show that the topological properties of the ELF in 2D are considerably 
simpler than in molecules studied previously. We compute the ELF and demonstrate its usefulness 
for various physical 2D systems, focusing on semiconductor quantum dots that effectively correspond 
to a confined 2D electron gas. The ELF visualizes the shell structure of harmonic quantum dots 
and provides insight into electron bonding in quantum-dot molecules. In external magnetic fields 
the ELF is found to be a useful measure of vorticity when analyzing the properties of quantum-Hall 
droplets. We show that the current-dependent term in the ELF expression is important in magnetic 
fields. 



I. INTRODUCTION 

The study of electron localization is linked to the pur- 
suit of a rigorous systematization of the elusive concept 
of a chemical bond^^ ubiquitous in quantum chemistry 
ever since the first attempts of its description^ The main 
point in the concept of electron localization is not where 
the electrons are (or in quantum terms where the elec- 
trons are likely to be), which can be monitored from the 
electron density, but where the electrons are localized. 

The concept of localization (and delocalization) arises 
from Pauli's exclusion principle: two electrons with the 
same spin cannot occupy the same spatial position. The 
consequence is the appearance of the Fermi hole^ a func- 
tion which determines how the probability of finding an 
electron at some point diminishes because of the pres- 
ence of another like-spin electron in its vicinity. In other 
words, the Fermi hole is that part of the exchange and 
correlation hole that stems solely from the fermionic na- 
ture of the electron and not because of the Coulomb re- 
pulsion. Any electron at a given point carries its asso- 
ciated Fermi hole as a counterpart. The localization or 
delocalization of the electron is equivalent to the localiza- 
tion or delocalization of its Fermi hole.^'^ Naturally, a lo- 
calization descriptor should be based on this function or, 
correspondingly, on the like-spin pair probability func- 
tion - this will be clarified below. Indeed, the so-called 
electron localization function (ELF)ii^i^ii^ is nothing else 
than an appropriate renormalization of the Fermi hole 
curvature 

The ELF is large where electrons pair. Our intuition 
in chemistry tells us that electrons should pair forming 
localized groups in each atomic shell of the inert cores, in 
the chemical bonds, and in the non-bonding or lone pairs. 
More precisely, the topological properties of a good local- 
ization function should partition the space for each group 
of electrons - also for delocalized groups of electrons such 
as TT systems in conjugated molecules, for which the ELF 
values should be low. Our understanding of chemistry 
is founded on the electron pair^ a consequence of Pauli's 
exclusion principle, and more particularly of the local- 
ization of one electron of each spin in a region of space. 



Accordingly, it is not surprising that the main purpose of 
the analysis of the ELF has been to help our intuition of 
the fundamental chemistry concepts of pairs and bonds. 

Due to this focus on the elucidation of the chemical 
bond, all studies of electron localization functions per- 
formed to date have referred to three-dimensional (3D) 
systems. To our knowledge, no attempts have been made 
to map the electron localization in two dimensions (2D). 
However, significant advances in semiconductor technol- 
ogy have enabled the production and manipulation of 
low-dimensional nanoscale structures. Common exam- 
ples of these systems are quantum dots (QDs)^ often 
also called artificial atoms. The electrons in QDs are 
confined in the interface of a semiconductor heterostruc- 
ture, e.g., GaAs/AlGaAs, so that the transverse dimen- 
sions controlled by a lateral confinement are considerably 
larger than the thickness of the QD. In most cases a 2D 
model describes the movement of the electrons with a 
reasonable accuracy. On the other hand, the approxi- 
mation for the in-plane confinement depends on the QD 
shape that can be, e.g., circular, rectangular, ring-like 
(quantum ring), or consist of several potential minima 
(quantum-dot molecule). Due to the tunable shape, size, 
and number of electrons, QDs have emerging applications 
in fields of quantum computation, and from a theoreti- 
cal point of view they are an ideal playground to study 
many-electron phenomena and test computational meth- 
ods. Hence, we expect the ELF to be of interest for 
various QD systems, since it carries valuable information 
about the electronic structure regardless of the dimen- 
sionality. 

This article is organized as follows. In Sec. [Til we 
present a detailed derivation of the ELF in 2D, focus- 
ing on the definition of the Fermi-hole curvature which 
is the core of the ELF, and show how the relevant ex- 
pressions hold compared with the 3D case. We also 
discuss the topological properties of the ELF in 2D. In 
Sec, mil we briefly present the QD model and our compu- 
tational scheme based on spin-density-functional theory. 
In Sec. IIVI we provide four examples where the ELF is 
a useful tool when analyzing the electronic properties of 
QDs. They include (i) visualization of the shell structure 
of a parabolic (harmonic) QD, (ii) showing the bond-like 
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interdot electron couplings in QD molecules, (iii) intro- 
ducing the ELF as a measure of vorticity in quantum-Hall 
droplets, and (iv) visualization of vortex localization in 
high magnetic fields. We also show that the current- 
dependent term in the ELF (or in the Fermi-hole cur- 
vature), which has been neglected in most studies until 
now, is relevant in order to obtain meaningful results for 
QDs in magnetic fields. Concluding remarks are given in 
SecEl 



II. ELECTRON LOCALIZATION FUNCTION 
A. Derivation of the ELF in two dimensions 

In the following we go in detail through the mathemat- 
ical derivation of the ELF. We focus on the differences in 
the expressions in 3D and 2D, in particular between the 
spherical and circular averages and kinetic-energy densi- 
ties. 

We consider the many-electron wave function 
^(ri<Ji, . . . ,rArcrAr), where N is the number of electrons, 
and r and a are the electronic position and spin coordi- 
nates, respectively. We point out that ^ can be either 
an eigenstate of a static Hamiltonian, or an evolving 
time-dependent state. Nevertheless, we omit the time 
coordinate in the notation. The necessary ingredients 
in the derivation of the ELF are the first-order reduced 
density matri x -'^^i-'^^i-'^^ 

Ci(ri|r'i) = 7V Yl /dr2...drA,x 

^*(ricri, . . . ,rArcrAr)^(ricri, . . ..Tnctn) , (1) 

and the diagonal of the second-order reduced density 
matrix^ 

N{N-1) ^ I dra . . . dr AT x 



(2) 



The diagonal of F^^^ is commonly known as the spin- 
density. 



^(r)=r(ij(r|r). 



(3) 



and the total density, n(r) = Xlcr^cr(^)' sums over the 
spin variable. It has an easy interpretation: n(r)dr is 
the probability (normalized to N) of finding an electron 
in a small volume dr at position r. The pair probability 
density is defined analogously by 



.(2) 



(ri,r2) = 7^^^)^^(ri,r2) 



(4) 



and can be interpreted in the following way: 
n*^^^(ri, r2)dridr2 is the probability of finding one 



electron at ri and, simultaneously^ another electron at 

1-2. 

Were the electrons fully independent, this latter prob- 
ability would just be the product of the electron densi- 
ties: n*^^^(ri,r2) = n(ri)n(r2). However, their exchange 
(fermionic character) and Coulomb interaction reduce 
the pair probability by a value known as the exchange 
and correlation hole: 



(ri,r2) = n(ri)n(r2) + /i(ri,r2) . 



(5) 



The Fermi hole is the part of /i(ri,r2) which is entirely 
due to the antisymmetric character of the wave function, 
regardless of the Coulomb interaction. Accordingly, it 
only appears for like-spin electrons. It is the dominant 
part at short distances (r2 ^ ri). 

The key function to study the electron localization is 
the like-spin conditional pair probability function defined 



as 



Pa(ri,r2) 



7ii^(ri,r2) 



(6) 



The physical meaning of this function is the following. 
It gives the probability of finding one a-spin electron at 
r2 knowing with certainty that another a-spin electron 
is at ri. The question is now how probable it is to find 
one like-spin electron in the vicinity of the first reference 
one. For this purpose it is useful to define a conditional 
pair probability function at a distance s via a spherical 
(in 3D) or circular (in 2D) average. Here we deal with 
the latter case, 



1 

p^(ri,5) = — / d(9 P^(ri,ri + su^) , (7) 
^TT Jo 

where uq = cos 6^x1 + sin 6^x2 . Since we are particularly 
interested in the behavior of Pcr(r, s) at small s, we need 
to perform Taylor expansions, 

^(ricr, (ri + sue)a, . . . , tnctn) = 

s {ue ' V2)^(ricr, r2cr, . . . , rArCrAr)r2=ri + 0{s'^) , (8) 

where V2 is the gradient with respect to the second elec- 
tron variable in ^. The term is absent due to Pauli's 
exclusion principle. We now take the square, 

|^(ricr, (ri + sue)a, . . . .tnctn)]'^ = 

I (u^ • V2)^(ria, r2a, . . . , rNCTN)r,=r^ f + 0{s^) = 
2 

J2 <^M(^>ii(^) + 0{s') , (9) 



where 



Ci = (Xi • V2)^(ricr, r2cr, . . . , rArcrAr)r 



(10) 
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with i = 1,2. The circular average is 

27T 



1 

2^ Jo 



dO |^(ricr, (ri + sue)(T, . . . , tnctn) 



2 -i p27r 



dO Ui{e)uj{0) = 



2 ^ 

i=l 



-S^ Vi|*(riC7, r2Cr, . . . , rN(JN)\l^=r, 



Ois') (11) 



This expression may then be used in combination with 
Eqs. (0), ®, and fT]) to obtain 



a(ri) = 



1 Vi7i'j(ri,r2). 

2 n,(ri) 



(12) 
(13) 



These equations are similar to the ones obtained in the 
3D case, in fact, they are the same except for the factor 
1/2 in Eq. (O, which is 1/3 in 3D. The function C^(r) 
satisfies exactly the same Eq. (fT3]) in 2D and in 3D. This 
is the function that measures the local like-spin pair prob- 
ability, is an inverse measure of electron localization, and 
is the function used to define the ELF. We should note 
also that V27aa (ri, r2)r2=ri is the Fermi hole curvature 
as defined by Dobson.^^ Other definitions, equivalent re- 
garding the information that they contain, are also pos- 
sible^ 

We conclude that the expressions that define the ELF 
do not change significantly from 3D to 2D. The func- 
tion Ccr(r) alone does not reveal the localization explic- 
itly. For that purpose one needs to perform the following 
renormalization, which defines the ELF as^ 



ELF,(r) 



1 



l + {a(r)/CHEG[^^(r)]} 



2 ' 



(14) 



where C^^^[ncr(r)] is the value of the Co- function for the 
homogeneous electron gas (HEG) of (constant) spin den- 
sities ncr(r). Analogously to the 3D case, the expression 
for C^^^ in 2D is the kinetic-energy density. 



HEG 



rio 



27rna 



(15) 



In terms of the total density n and the polarization (often 
called also magnetization) ( = (n^ — ni)/n, the expres- 
sion immediately yields 



c."'^°KC] = ^(i + C> 



(16) 



for the total kinetic-energy density. 

The above definition of Co- is completely general, but 
requires the diagonal of the second-order reduced density 
matrix, which is a fairly complex and often unmanageable 
object. In most cases, however, a further approximation 



is taken by assuming that ^ is a single-determinantal 
wave function. In practice, this means that one uses 
the Hartree-Fock approximation, therefore neglecting the 
effect of correlations in the ELF , or one uses density- 
functional theory (DFT) (see below), and calculates the 
ELF of the non-interacting auxiliary Kohn-Sham (KS) 
system. Both possibilities are in principle unsatisfactory, 
but the experience has demonstrated that the results are 
usually unaffected despite a notable mathematical sim- 
plification (see Ref. [ly] for the detailed steps). If ^ is 

formed by the set of KS orbitals {v^i^j^Ji and 

{N^ + A/"; = TV) within the spin-DFT (SDFT), one finds 



a(r)=r,(r) 



l (Vn.(r))^ 
4 n^(r) n^(r) 



(17) 



where ncr(r) = Ylf=i |v^icr(r)p is the spin-density [cf. 
Eq. ([3])], To- is the kinetic-energy density given by 



r,(r) = ^|V(p,,(r)|^ 



(18) 



and jp^cr is the spin-resolved paramagnetic current density 

ipA^) = ^Y1 [(^^- W) W - ^*a(r)Vv:)..(r)] . 

(19) 

Here we point out that most ELF calculations refer to 
closed-shell molecules (or finite systems in general) in 
their ground state and in the absence of magnetic fields 
or the spin-orbit coupling. In those cases real wave func- 
tions can be assumed and the current term in Eq. (pT|) is 
absent. In Sees. IIV Cl and llV Dl we note the importance of 
the current term, which is needed in the ELF also during 
time-dependent processes^^i^ 



B. Topological properties 

The basis for the topological investigation of the ELF 
was established by Silvi and Savin. ^ The information con- 
tained in the ELF is extracted with the help of some 
basic concepts borrowed from the theory of dynamical 
systems.^!^ The ELF is a scalar continuous function 
bounded between and 1, and one should look at its 
gradient field, which in turn defines the set of attractors 
- roughly speaking the maxima of the ELF. The impor- 
tant concept is then that of a basin of each attractor - 
the set of points for which that attractor is the cj-limit. 
In chemistry, each basin is then identified either as a core 
basin (it contains a nucleus), or as a valence basin. The 
latter ones are then classified according to the number of 
core basins that they have frontier with: if they have only 
one neighbor core, they are called monosynaptic; if they 
have more than one (disynaptic, trisynaptic etc.), they 
are bonds. The absence of real nuclei in the 2D models 
makes this distinction not relevant. However, the divi- 
sion of space into localization basins is still pertinent: the 
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basins give us information about the groups of electrons; 
exchange of electrons is more unlikely between different 
basins than inside them. In 3D the attractors can then 
be single points in the absence of special symmetries, 
ring-shaped lines of attractors around a symmetry axis, 
or spheres around isolated atoms. In 2D, obviously, the 
topology simplifies and we can only find, besides point 
attractors, ring attractors for cylindrically symmetrical 
problems. 



III. QUANTUM-DOT MODEL 



We study the ELF in 2D QDs restricted to the xy plane. 
The system is described by an effective-mass TV-electron 
Hamiltonian 



H 



N 



2m* 



Ve^tiVi) ^ Ez,i 



Kj 



(20) 

where A is the external vector potential (in symmet- 
ric gauge) of the homogeneous magnetic field 'B = Bqz 
perpendicular to the plane, lext(r) is the external con- 
fining potential in the xy plane (see below), and Ez = 
9* I^bSzBo is the Zeeman energy. We apply the conven- 
tional effective- mass approximation with the material pa- 
rameters for GaAs: the effective mass m* = 0.067, the 
dielectric constant e* = 12.4, and the gyromagnetic ratio 
^* = -0.44. 

We solve the ground-state problem associated with the 
A/'-electron Hamiltonian ([2Q|) by applying the SDFT in 
the collinear-spin representation. The KS states, needed 
in Eqs. (fT7|) - (p!9|) for computing the ELF as defined in 
Eq. (dH), are solved from the KS equation 



(-zV. + ^A(r))2 



2m* 



+ ^Ks(r) 



ipia{v) = eia(pia{r), (21) 



where the KS potential V-^^{r) is a sum of the ex- 
ternal confining potential, the Hartree potential, and 
the exchange-correlation potential given by V^^(r) = 
(5£^xc[^^7 ^^]/^^^(r). To approximate the exchange- 
correlation energy E^c we use the local spin-density ap- 
proximation (LSDA) with a parametrization provided by 
Attaccalite and co-workers. In QD systems the SDFT 
scheme together with the LSDA leads to good numeri- 
cal accuracy in comparison with quantum Monte Carlo 
calculations, even in relatively high magnetic fields 
We also point out that the local- vorticity approximation 
within the current-SDFT^^ does not lead to a consider- 
able improvement over the SDFT results.^^ In the numer- 
ical calculations we apply the real-space octopus code^^^ 
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FIG. 1: (color online). ELF (red solid lines) with the cor- 
responding electron densities (blue dashed lines) for two- 
dimensional closed-shell quantum dots. The maximum den- 
sity values are scaled to one. 



IV. EXAMPLES 
A. Shell structure 

First we examine the ELF in parabolic QDs at zero 
magnetic field {B = 0). Thus, we choose the exter- 
nal confining potential in Eq. ([2Q|) as Ve^t{r) = uJ^r^ 
where we set the confinement strength to uj^ = 0.42168 
a.u. =5 meV. The system is a 2D harmonic oscillator 
with single-electron eigenenergies Cni = (2n + |/| + l)u;o^ 
where n = 0, 1, 2, . . . is the radial and / = 0, ±1, ±2, . . . 
is the azimuthal quantum number. The corresponding 
shell structure has been experimentally well depicted in 
the addition energies of a double-barrier GaAs QD^ 
Further details in the measured addition energies result- 
ing from the electron-electron interaction, such as Hund's 
rule -type of behavior, have been theoretically verified in 
numerous studies employing, e.g., SDFT^^i^i^ 

Figure [T] shows the ELF (red solid lines) and the cor- 
responding electron densities (blue dashed lines) when 
TV = 6, 12, 20, 30, 42, and 56. In QDs these "magic" elec- 
tron numbers correspond to the first closed-shell config- 
urations when N > 2. The shell structure is visualized 
clearly by the ELF. We note, however, that the local max- 
ima in the ELF do not directly correspond to the filled 
energy shells. Instead, each peak can be associated with 
a doubly occupied single-electron state on the highest en- 
ergy shell. This means that the probability of electron 
localization is highest close to the Fermi level. For exam- 
ple, in a 12-electron QD shown in Fig. [TJb), the highest 
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energy shell has six electrons with / = 0, ±2, leading to 
distinctive peaks in the ELF at the center (corresponding 
to / = 0) and at r ~ 4 (corresponding to / = ±2). Nev- 
ertheless, in this system the number of energy shells is 
equal to the number of extrema in the ELF. As expected, 
the variation between the extrema decreases as a function 
of A^, since the highest shell becomes relatively less dom- 
inant. As seen in Fig.[Tl the shell structure is visible also 
in the electron density n(r), but it is less clear than in 
the ELF, especially at large electron numbers where the 
structure is barely visible. The ELF instead provides an 
unambiguous visualization of the shell structure, which is 
well plausible from its physical nature discussed in Sees. [J 
andini 



B. Quantum- dot molecules 



The applicability of the ELF to visualize pairs and 
bonds in molecules immediately suggests to use the 
ELF for the 2D counterparts commonly known as QD 
molecules (QDMs). Since the spin-qubit proposal of Loss 
and DeVinzenzo^ in 1998, coupled QD systems have at- 
tracted wide interest both experimentall}*^ and theoret- 
icall}*^, particularly in terms of charge and spin manipu- 
lation. Following the standard QDM definition, we write 
the external confining potential as 



Ve^t{r) -uJq min 



M 



(22) 



where M is the number of potential minima located at 
Tj = (xj^i/j). The case of M = 1, = (0,0) is equal 
to a single harmonic QD considered above. Here we set 
M = 4 and Vj = (±V^, ±V2) with = 0.5 a.u. This 
corresponds to a square-symmetric QDM with four min- 
ima^^ Now, Fig. E] shows the result of our SDFT cal- 
culations for the total densities and ELFs when A" = 8 
(a-b). A" = 12 (c-d), and A" = 16 (e-f), respectively. The 
ground-state total spin 6* = in all the cases shown. The 
densities of A^ = 8 and A' = 12 QDMs look very similar 
with four distinctive peaks, whereas the A" = 16 QDM 
shows more structure around the density maxima. How- 
ever, the ELFs of these QDMs are very different from 
each other, having four maxima with one (local) mini- 
mum (A^ — 8), five maxima with four minima (A^ = 12), 
and four maxima with five minima (A^ = 16). Despite 
the fact that the concept of chemical bonding seems less 
meaningful in systems where the confining potential is 
fixed (such as in 2D QD-systems), the central "basin" 
(see Sec. IHB[) in Fig. [2fd), for example, is very infor- 
mative: there is pronounced localization at the center, 
where the confining potential actually has a repulsive 
cusp. Furthermore, this ground-state property is not at 
all visible in the bare electron density shown in Fig.[2](c). 
Therefore, we find that the ELF reveals features in the 
electronic structure of QDMs which are absent in the 
electron density - a fact also demonstrated in 3D. 



FIG. 2: (color online) Ground-state total densities and ELFs 
for four- minima quantum- dot molecules having 8, 12, and 16 
electrons, respectively. The total spin 5* = in all the cases. 



C. Maximum-density droplet 

Next we perform our analysis of the 2D ELF for non- 
zero external magnetic fields {B > 0). As the first exam- 
ple we consider a maximum-density droplet (MDD) ^^i^^ 
of a single parabolic QD. The MDD is a fully-polarized 
{S = N/2) state, where the electrons have consecutive 
angular momenta from / = 0to/ = — A' + l. Each elec- 
tron is accompanied by a so-called Pauli vortex which 
corresponds to a change of 27r in the phase of the many- 
electron wave function. Hence, the MDD state is a unique 
finite-size counterpart of the filling-factor u = 1 quantum 
Hall state of the uniform 2D electron gas. Figures [3](a) 
and (b) show the electron density and the ELF of the 
MDD state in a 12-electron parabolic QD at B = 8 
T. The MDD density has the well-known fiat shape^i^ 
whereas the ELF is characterized by a fiat interior and 
localization around the edge of the QD. Intersections of 
the density and the ELF are plotted in Figs. [3](c-d) (red 
dashed lines), together with the results for 6 and 20- 
electron QDs at 5 = 6 T and 8 T, respectively. Inter- 
estingly, the values of the ELF are very close to 1/2 in 
the interior of the MDDs. Thus, from Eq. ([T4j) we find 
Ccr(r) ^ C™^(r) in this regime, i.e., the kinetic-energy 
density of the uniform 2D electron gas. The result sug- 
gests that the localization of electrons and vortices com- 
pensate each other, which is well plausible considering 
the V = 1 character of the MDD state. We may thus 
expect that at higher magnetic fields, corresponding to 
the fractional quantum-Hall regime {v < 1), ELF values 
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FIG. 3: (color online) (a,b) Electron density and ELF of the 
maximum- density droplet in a 12-electron quantum dot. (c,d) 
Intersections of (a,b) together with the results for = 6 and 
= 20 quantum dots. (e,f) Paramagnetic current densi- 
ties and the kinetic-energy densities of the maximum-density- 
droplet states, respectively. 



below 1/2 can be found at positions of high vorticity. 

Figure [S^d) shows also the ELF of the N = 20 QD cal- 
culated without the current-density term jp in Eq. (pT|) . 
That curve equals to the full ELF at r = where the cur- 
rent density is zero, but decays exponentially at larger 
r. Obviously, that result does not capture the correct 
behavior of the MDD state. The dramatic difference 
demonstrates the importance of the current-density de- 
pendence in the ELF expression, already at relatively 
moderate magnetic fields. In Fig. [Sje) we plot jp which 
increases linearly due to the successive increase in the an- 
gular momenta of the KS states in the MDD (see above). 



D. Vortex structures 

Increasing the magnetic field above the MDD limit at 
= 1, which, as seen above, corresponds to ELF=l/2, 
leads to localization of vortices. Depending on the QD ge- 
ometry, the vortices may form clusters^Si^i^i^ or merge 
together yielding multiple phase quantization!^ In Fig.[4l 
we show the electron density and the ELF of a single- 
vortex solution in a six-electron parabolic QD at 5 = 11 
T. In this case the vortex is localized at the center and 
directly visible as a hole in the density. As expected, the 
ELF shows a similar strucure. The ELF has a bump 
with a value ~ 1/2, located between the center and 
the edge of the QD. This region can be interpreted to 
have, on the average, a local balance between electrons 
and (Pauli-)vortices, and it separates the localized vortex 
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FIG. 4: (color online). Electron density (a) and the ELF (b) 
of a single- vortex solution in a six-electron parabolic quantum 
dot at 5 = 11 T. The lower panel shows middle intersections 
of the figures in the upper panel. 




X (a.u.) X (a.u.) 



FIG. 5: (color online). Electron density (a) and the ELF (b) 
of a two- vortex solution in a six-electron rectangular quantum 
dot at 5 = 16 T. The lower panel shows intersections of the 
upper figures at yi and ^2, marked in the upper right figure. 



(ELF= 0) from the edge having high electron localization 
(ELF^ 0.8). 

Increasing the magnetic field further leads to the for- 
mation of more vortices. Figure [5] shows a two- vortex 
structure in a six-electron rectangular (hard-wall) QD 
at 5 = 16 T. The system has been studied in detail in 
Refs. [Ill and [3]. Again, the vortices are seen as zeros in 
the ELF, whereas the density is not exactly zero at these 
positions. Actually, the numerically exact density at the 
vortices is even further above zero than the SDFT den- 
sity due to the multiconfigurational nature of the many- 
electron wave function. Nevertheless, the ELF=0 result 
at the vortex positions is plausible, reflecting again the 
difference between the probable location (density) and 
the localization (ELF) of the electrons as discussed in 
Sec. Ill Along the edge of the QD, the ELF in Fig. E^b) 
shows six clear peaks with maximum values close to one, 
separated by local minima where ELF~ 1/2. Thus, in 
this QD the ELF reveals the Wigner crystallization, i.e., 
localization of electrons around their classical positions 
which are determined by the geometry of the system. 
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This effect is considerably less pronounced in the elec- 
tron density shown in Fig. [51(a). 

V. CONCLUDING REMARKS 

The electronic structure of a many-electron system is 
fully characterized by its many-body wave function. In 
order to acquire an intuitive visual understanding of the 
system, however, we must look at simpler objects - inte- 
grated magnitudes such as the electronic density, which 
lives in a lower dimensional space. Unfortunately, the 
density does not fully reveal all the intricacies of the elec- 
tronic structure, even if we know that it contains all the 
information. 

In the past, the ELF has proven to be a useful com- 
panion to the density in the task of providing us with in- 
sightful intuition on the electronic structure of molecules. 
In this work we have defined the ELF in 2D, and we have 
demonstrated that its visualization, in addition to that 
of the density, helps to understand the electronic struc- 
ture of 2D systems such as semiconductor quantum dots. 
We have shown the usefulness of the ELF to visualize 
the shell structure, as well as the bond-like features in 



coupled systems. In particular, we have found that in 
magnetic fields the ELF can be used as a measure of vor- 
ticity, revealing the local relation between the localization 
of electrons and vortices in the system. In this context, 
we have shown that the current-dependent term, which 
has been neglected in previous 3D studies, has a major 
contribution to the ELF in magnetic fields. We expect 
that alongside the rapid technological developments in 
the fabrication and manipulation of low-dimensional sys- 
tems, the ELF will prove to be a universally applicable 
theoretical tool to obtain detailed information of both 
static and dynamic many-particle properties. 
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